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We give a physical picture of the low-energy sector of the spin-i kagome Heisenberg antiferro- 
magnets (KAFs). It is shown that kagome lattice can be represented as a set of blocks containing 
12 spins, having form of stars and arranged in a triangular lattice. Each of these stars has two 
degenerate singlet ground states which can be considered in terms of pseudospin-i. It is shown 
using symmetry consideration that the KAF lower singlet band is made by inter-star interaction 
from these degenerate states. We demonstrate that this band is described by effective Hamiltonian 
of a magnet in the external magnetic field. The general form of this Hamiltonian is established. The 
Hamiltonian parameters are calculated up to the third order of perturbation theory. The ground 
state energy calculated in the considered model is lower than that evaluated numerically in the 
previous finite clusters studies. A way of experimental verification of this picture using neutron 
scattering is discussed. It is shown that the approach presented can not be expanded directly on 
KAFs with larger spin values. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Gb 



I. INTRODUCTION 

Unusual low-temperature properties of kagome anti- 
ferromagnets (KAFs) attracted much attention of both 
theorists and experimenters in last decade. Seemingly 
the most intriguing features were observed in specific heat 
C measurements of SrCrgpGai2-9pOi9 (spin-| kagome 
material) jj- It has been revealed a peak at T sa 5 K 
which was practically independent of magnetic field up to 
12 T and C appeared to be proportional to at T < 5 
K. 

There was no appropriate theory describing KAF low- 
energy sector so far. Qualitative understanding of the 
KAF low-temperature physics was based mostly on re- 
sults of numerous finite cluster investigations |^ H, ^ |M 
1^ . They revealed a gap separating the ground state from 
the upper triplet levels and a band of nonmagnetic sin- 
glet excitations with a very small or zero gap inside the 
spin gap. The number of states in the band increases 
with the number of sites N as a^. It was obtained for 
samples with up to 36 sites that a = 1.15 and 1.18 for 
even and odd N, respectively It is believed now 

that this wealth of singlets is responsible for specific heat 
low-T peak and explains its field independence P, 01 • 

The origin of the singlet band as well as the nature of 
the ground state were unclear until now. Previous exact 
diagonalization studies of clusters with A^ < 36 0, Q 
revealed exponential decay of the spin-spin and dimer- 
dimer correlation functions. So the point of view that 
KAF is a spin liquid is widely accepted 011 H H 0, II 

UMMm 

It seems the best candidate for description of KAF low- 
energy properties is a quantum dimer model 6^ J^J . In 
Ref. [iJI an approach firstly pioneered by Subrahmanyam 
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in Ref. |15| was developed in which a spin-i kagome lat- 
tice is considered as a set of interactive triangles with a 
spin in each apex. It was suggested there to work in the 
subspace where the total spin of each triangle is 1 /2 inves- 
tigating low-lying excitations. In spite of certain success 
of this approach (coincidence of the low-energy spectrum 
and the number of lower singlet excitations in samples 
with up to 36 sites with exact diagonalization results) a 
further development of this approach is required to get 
the full physical description of KAF. 

Recently in Ref. it was suggested a model of frus- 
trated antiferromagnet which low-energy properties can 
be generic for KAF as well as for two other types of 
frustrated magnets which possess a similar behavior as 
KAF and have many singlet states inside the triplet gap: 
pyrochlore ^] and CaV40g 0|. Weakly interactive 
plaquettes on the square lattice were considered there. 
Each plaquette has two almost degenerate singlet ground 
states, so a band of singlet excitations arises if the inter- 
plaquette interaction is taken into account. It is shown 
that there is a quantum phase transition in the model at 
a critical value of frustration separating a disorder pla- 
quette phase and a columnar dimer one. In the proximity 
of this transition the specific heat has a low-T peak below 
which it possesses a power low-temperature dependence. 

In our recent paper Ref. we have shown that a sim- 
ilar picture is relevant for spin-i KAF. It was proposed to 
consider a kagome lattice as a set of stars with 12 spins 
arranged in a triangular lattice (see Fig. Numerical 
diagonalization has shown that star has two degenerate 
singlet ground states separated from the upper triplet 
levels by a gap. These states form a singlet energy band 
as a result of inter-star interaction. It was assumed that 
this band determines the KAF low-energy singlet sector. 
We have shown that it is described by Hamiltonian of 
a magnet in the external magnetic field where degener- 
ate states of the stars are represented in terms of two 
projections of pseudospin-i. 
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The main advantage of a star over a triangle proposed 
m Refs. MM 

as a starting point to study KAF is that 
consideration of star two singlet ground states explic- 
itly leads to effective Hamiltonian which describes singlet 
low-temperature dynamics. 

In the present more comprehensive paper we develop 
this approach. It is proved using symmetry considera- 
tions presented in Sec. ^1 that the singlet band above 
discussed determines the KAF lower singlet sector. 

This band is studied in Sec. IIIII The general form 
of the effective Hamiltonian is established there. The 
Hamiltonian parameters are calculated up to the third 
order of perturbation theory. The ground state energy 
calculated in the considered model is lower than that 
evaluated numerically in the previous finite clusters stud- 
ies. We demonstrate that our approach can not be ex- 
panded directly on KAFs with larger spin values. A way 
of experimental verification of this picture using neutron 
scattering is discussed. 

Finally we summarize our results in Sec. IIVI 



II. SYMMETRY CONSIDERATION 

We start with the Hamiltonian of the spin-i kagome 
Heisenberg antiferromagnet: 



functions can be represented as follows: 

1 



Ti-o — Ji S^Sj + J2 SiSj 



(1) 



where and denote nearest and next-nearest 

neighbors on the kagome lattice, respectively, shown in 
Fig. n The case of IJ2I ^ -/i is considered in this paper. 
We discuss a possibility of both signs of next-nearest- 
neighbor interactions — a ferromagnetic and an antifer- 
romagnetic one. As is shown below, in spite of the small- 
ness the second term in Eq. iQJ can be of importance for 
the low-energy properties. 

Kagome lattice can be represented as a set of stars 
arranged in a triangular lattice (see Fig. To begin 
with we neglect interaction between stars and put J2 = 
in Eq. A star is a system of 12 spins. Let us consider 
its properties in detail. 

As the Hamiltonian Eq. Q commutates with all the 
projections of the total spin operator, all the star levels 
are classified by the values S, irreducible representations 
(IRs) of its symmetry group and are degenerated with 
Sz. The star symmetry group Cg^ contains six rotations 
and reflections with respect to six lines passing through 
the center. There are four one-dimensional and two two- 
dimensional IRs which are presented in Appendix ^ In 
their basis the matrix of the Hamiltonian has a block 
structure. Each block has been diagonalized numerically. 

As a result it was obtained that the star has doubly de- 
generated singlet ground state separated from the lower 
triplet level by a gap A « 0.26Ji. Ground state wave 



v/2 + 1/16 
1 

V2-1/16 



(01 - 02), 



(2) 
(3) 



where functions 0i and 02 sue shown schematically in 
Fig. 13 A bold line there represents the singlet state of 
the corresponding two spins, i.e. (|T)Ji)j - |i)ilT>j)/V2- 

It can be shown that 0i and 02 are not orthogonal: 
their scalar product is (0i02) — 1/32. They contain six 
singlets each having the energy ~S{S + l)Ji = — 3/4Ji. 
One can obtain that interaction between singlets does 
not contribute to the energy of the ground states which 
is consequently equal to — 4.5Ji. 

Functions 0i and 02 are invariant with respect to ro- 
tations of the star and they transform into each other 
under reflections. Hence is invariant under all the 
symmetry group transformations. In contrast, function 
is invariant under rotations, changes the sign with 
reflections and is transformed under the representation 
(|A3I) . So the ground state has accidental degeneracy. As 
it is shown in the next section, the next-nearest neighbor 
interaction, which has the same symmetry as the original 
Hamiltonian, lifts this degeneracy. 

KAF containing M non-interactive stars has an energy 
spectrum with great levels degeneracy when A/" 3> 1 . For 
example, the ground state degeneracy is 2^ and that of 
the lowest triplet level is 3A/'2*^^~^'. Interaction between 
stars gives rise to an energy band from every such group 
of levels and it is a very difficult task to follow their 
evolution. Meanwhile group theory allows one to make 
some conclusions about the KAF low-energy sector. We 
show now that the singlet band stemmed from the ground 
state can not be overlapped by those originated from the 
upper singlet levels. 

Let us consider a cluster with seven stars shown in 
the Fig. neglecting for a beginning interaction be- 
tween them. The symmetry group of the cluster is Cqv 
also. The ground state has degeneracy 2^ — 128. Cor- 
responding wave functions transformed under IRs of Cgu 
are constructed as linear combinations of products con- 
taining or for each star. It is easy to show using 
a standard procedure of irreducible representations ba- 
sis construction [23, |^ and Appendix ^ that there is 
at least two ground state wave functions of the cluster 
discussed transformed under any given IR. 

It is important to mention that operator of inter-star 
interaction in the cluster has the same symmetry as the 
intra-star one which is a sum of the stars Hamiltonians. 
Then the interaction between stars commutates with the 
square of the total spin operator. So if we begin to in- 
crease the inter-star interaction from zero, all the levels 
would move by energy but their classification can not be 
changed. Levels can cross with each other as the interac- 
tion arises from to Ji but the crossing is forbidden for 
those of the same symmetry. It is a consequence of the 
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symmetry theorem proved in Refs. |20l |21|. Hence, one 
can conclude that the lower singlet sector of the cluster 
is formed by states which stem from the origin 128 lower 
levels. 

We can lead to the same conclusions considering clus- 
ters of the symmetry Cgu with the greater number of 
stars. So we will assume in the following that the KAF 
low-energy singlet sector is formed by the states origi- 
nated from those in which each star is in one of the states 
^f^- or 4"-. Because bands with S ^ Q can overlap the 
singlet one under discussion we have to suppose that the 
KAF low-energy properties are determined by the lowest 
states in this singlet band. 

As interaction between stars commutates with the 
square of the total spin operator, bands of the different 
iS" can be studied independently. The KAF lower singlet 
sector is considered in detail in the next section. Inves- 
tigation of states with 7^ is out of the scope of this 
paper. 



III. SINGLET DYNAMICS 

In this section we derive the general form of the ef- 
fective Hamiltonian described the lower singlet sector. 
Inter-star interaction is considered as a perturbation. De- 
spite of its not being small comparing to the intra-star 
one there are reasons presented below to use perturbation 
expansion here. 

Two-stars coupling. — We start with consideration of 
the interaction between two nearest stars still neglecting 
the second term in Eq. Initially there is four times 

degenerated ground state with wave functions {'^n}'^n}} 
(here = — and the upper index labels the stars) and 



the energy En}n2 



(0) 



-9Ji. As it is seen 



from Fig. 01 the interaction has the form 



V = Ji{SY>S 



(1)«(2) 



c(i)c(2) 



(4) 



According to the standard theory [23 the followmg con- 
ditions should be fulfilled to consider y as a perturbation: 



I ^mim2 I 



k(0) _ 

^n\n2 ^m\m2 



AO) 



« 1, 



(5) 



,(2)n 



where Vn,n2;rmm2 = I V \ ^^^\ 

denotes excited singlet levels of the two stars and Ui ~ 
+,—■ We have calculated C^^^ for rii — using 
wave functions obtained numerically and found that all 
of these coefficients do not exceed 0.09. So conditions 
Eq. @ are fulfilled. Then the maximum value of the 
sum X^mims I^mim2l^ ^"^^ which is also small enough. 
Thus interaction between stars will be considered as a 
perturbation in the following. 

We proceed with calculations of corrections to the ini- 
tial ground state energy of two stars. For this purpose, as 



the state is fourfold degenerated, one has to solve a secu- 
lar equation p3 |. The corresponding matrix elements in 
the third order of perturbation theory are given by po| 



^n\n2\'m\ra2 ^raira2\kik2 



p(0) 



(0) 



E E 



mi,m2 gi,g2 



mi,m2 ^nin2 ^mim2 
^iiTi2 ■,mim2 "12:9192 ^qiq2\ki k2 



(0) 



AO) 



"2 



pj(0) ^ 

^9192 ) 



, (6) 



h, — . Obviously the first term in Eq. © 
is zero and the second one can be represented as follows: 



where rii , ki 



-I j dte~ 
'0 



Hnin2'-ki k2 

X (*«*(2) I ye-*(<'+^^")F I (7) 

where ^ are Hamiltonians of the corresponding stars. 
The third term in Eq. 10 will be considered later. Us- 
ing the symmetry of the functions and discussed 
above and invariance of the system under reflection with 
respect to the dotted line plotted in Fig. ^ one can show 
that only nonzero elements belong to the first and to the 
second diagonals (i.e. with ni = fci, n2 = fc2 and with 
TT-i 7^ fci, n2 7^ ^2). We have calculated them numerically 
with a very high precision by expansion of the opera- 
tor exponent up to the power 150. The results can be 
represented in the following form: 



H+-,+- 



-fli + 02 - as, 
-ai + as, 
-ai + 03, 
-fli - 02 - 03, 



(8a) 
(8b) 
(8c) 
(8d) 



where ai = 0.256Ji, 02 = O.OlSJi and 03 = 0.0017Ji. 

Terms of the second diagonal = H = 

—H^ ; y = —H |_;_| = 04 = — O.OOO2J1 are much 

lesser then ai, 02 and 03. So the interaction shifts all 
the levels on the value — ai and lifts their degeneracy. 
Constants 02 , 03 and 04 determine the levels splitting. It 
is seen that the splitting is very small comparing to the 
shift. 

So KAF appears to be a set of two-levels interacting 
systems and one can naturally represent the low-energy 
singlet sector of Hilbert space in terms of pseudospins: 
It) = and li) = It is seen from Eqs. © that in 
these terms the interaction between stars is described by 
the Hamiltonian of a ferromagnet in the external mag- 
netic field: 



c 



(9) 



where labels now nearest-neighbor pseudospins, ar- 
ranged in a triangular lattice formed by the stars, s is 
spin-i operator, C = -5.268Ji7V, J, = ia^ = -0.007 Ji, 



Jy = 4a4 



-O.OOlJi and h 



-602 



-0.092Ji. Here 
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M = N/12 is the number of stars in the lattice. Factor 
6 appears in the expression for h because each star in- 
teracts with 6 neighbors. One can see that the magnetic 
field in the effective Hamiltonian Eq. is much more 
larger than the exchange. So in this approximation stars 
behave like almost free spins in the external magnetic 
field and the ground state of KAF has long-range singlet 
order which settles on the triangular star lattice and is 
formed by stars in states. 

corrections. — The field remains the largest term of 
the effective Hamiltonian and the KAF ground state has 
the same long-range order if one takes into account 
terms in the perturbation series. For two-stars coupling 
corrections have the form given by Eq. ®. Terms 
imply analysis of the three-stars interaction as well. 
Nonzero contributions from them one gets only for the 
configuration presented in Fig.^Ib). The secular matrix 
for three stars is of the size 8x8. We have calculated 
corrections with a very high precision using the integral 
representation similar to that Eq. Q for the second term 
in Eq. © . All the operator exponents were expanded up 
to the power 150. As a result the low-energy properties 
of the KAF are described by the effective Hamiltonian: 

^ = E [^^'^'^1 + ^-'^'j + ^y^'"'] +hY,st+C, 

(10) 

where all parameters are given in the Tabled It describes 
two-stars coupling. We omit the three pseudospins terms 
in Eq. (|10|l which have the form s^Sjsl and sfs^s^ and 
describe three-stars interaction. The corresponding co- 
efficients are of the order of 10~^ Ji and lO""* Ji, respec- 
tively, and are negligible in comparison with those of the 
presented terms. It should be stressed that within our 
precision Hamiltonian Eq. (|10|l is an exact mapping of 
the original Hcisenberg model in the low-energy sector 
(excitation energy lj ~ max{ J^, J^, Jx,h} <^ Ji). 

As one can conclude from the Tabled studying cor- 
rections the common shift given by them remains much 
larger than the levels splitting in both cases of two- and 
three-stars coupling. Meanwhile the values of pertur- 
bation terms are approximately two times smaller than 
those of ones. So the change in the effective Hamil- 
tonian from terms is significant and analysis of the 
perturbation series can not be finished at this point for 
correct establishing of the effective Hamiltonian. Unfor- 
tunately such a work requires great computer capacity 
we do not dispose now. We have to restrict ourself with 
this precision here. 

One can judge about the applicability of the perturba- 
tion series from the following values of the ground state 
energy of two interacting stars shown in Fig. ^ obtained 
numerically and using first two orders of perturbation 
theory. Ground state energy of two non-interacting stars 
is — 9Ji. That of two interacting ones which we have cal- 
culated numerically by power method is — 9.62Ji. 
Whereas the ground state energy obtained using Ta- 
ble m is — 9.42 Ji (contribution from and terms 



are — 0.27Ji and — O.lSJi, respectively). 

Effective Hamiltonian structure. — Despite the pertur- 
bation theory works quite badly in the star model and 
many perturbation terms are to be taken into account, 
we can proof now that Eq. (|10|l is the most general form 
of the effective Hamiltonian assuming that n-pseudospins 
couplings with n > 2 are small as it was in the case of 
n = 3 discussed above. Let us consider possible terms 
of the form sfsf, sfs~, sf and s^ . In these cases num- 
bers of functions and to the right hand of the 
corresponding matrix elements differ from those to the 
left hand by unit. As it was also pointed out above, 
kagome lattice contains lines of symmetry reflections and 
star Hamiltonian and inter-star interaction are invariant 
with respect to these reflections. As 4'+ are invariant 
and 'J- changes the sign under these transformations, 
the matrix elements are equal to themselves with the op- 
posite sign and so must be zero. Other possible term 
sfsj = -(z/4)(s+s^ - s~s~ + s^sf — sfs~) can not 
appear in the effective Hamiltonian because the corre- 
sponding matrix elements should be imaginary. 

Ground state. — As is clear from the Tabled and h 
are the largest parameters of the Hamiltonian Eq. I|10|) 
in our approximation. So KAF behaves like Ising anti- 
ferromagnet in perpendicular magnetic fleld. In this case 
the classical value of the fleld at which spin-flip occurs 
is hs-f = Jx which is approximately 2.6 times smaller 
than h. So the ground state should remain ordered with 
all the stars in ^ state. 

The ground state energy and that of the upper edge of 
the singlet band calculated using Table n are (-4.5 Ji + 
AC + hjl + iJ^./A^N = -0.452 Ji TV and (-4.5 Ji AC - 
h/2 + iJz/A)!^ = -0.437JiiV, respectively. Corrections 
from Jx to these values in the flrst non-zero order of 
perturbation theory are given by {'i/\%)Nj^/h and are 
negligible. At the same time the ground state energy of 
the largest cluster with = 36 considered before nu- 
merically was — 0.438 JiA^ [3- So we beheve that clusters 
used in the previous studies were too small to reflect the 
Heisenberg KAF low-energy sector at J2 = properly. 

Interaction J2. — We show now that in spite of its 
smallness the next-nearest neighbor interaction can play 
an important role for low-energy properties. We have 
calculated J2 corrections to the parameters of the effec- 
tive Hamiltonian Eq. I|10|) for the first and for the second 
terms in Eq. ((HJ only. There are 12 intrinsic J2 inter- 
actions in each star which splits the doubly degenerated 
ground state and, as is seen from the Table |l| gives a 
contribution to the magnetic field h and to the constant 
C. 

As is clear from Fig. ^ the two-stars coupling is given 
now by the operator V = J2(S^^^S^^^ + S^^^sf^ + 

''Sg^-' + S3^''S2^'')- Corrections proportional to J2 were 
calculated by the same way as above and are presented 
m the Table |I| too. It is seen that the contribution of 
next-nearest interactions to the magnetic field becomes 
significant if | J2I ~ O.lJi. If J2 < (ferromagnet inter- 
action) they can even change the sign of h. 
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Effect of next-nearest ferromagnetic coupling for KAF 
properties was previously studied in Ref. 2] numerically 
on finite clusters with < 27 in the wide range of the 
values of J2. It was shown there that at \J2\IJ1 ~ 1 
the ground state has -s/S x VS magnetic structure. At 
1-^21/^1 <C 1 the ground state was found to be disorder 
and there is a band of singlet excitations inside the triplet 
gap. As we demonstrated above in our approach this 
band is a result of star ground state degeneracy. 

specific heat behavior. — There were many specula- 
tions on the low-temperature dependence of the KAF 
specific heat C (xT'^ observed experimentally for 5* = 3/2 
(see Refs. and references therein). As we obtained 

above, low-T properties are described by effective Hamil- 
tonian Eq. H10|l of a magnet which has spectrum of the 
form eq — \J [cqY + A'^ at g <C 1 and can be in ordered 
or disordered phases depending on particular values of 
the Hamiltonian parameters. Small A' here signifies the 
proximity to the quantum critical point in which C oc T^. 
Such a situation arises in the singlet dynamics of the 
model of interactive plaquets fl^ discussed above. We 
do not present here the corresponding analysis because 
parameters of the effective Hamiltonian could be changed 
in the further orders of the perturbation theory. 

Experimental verification. — In both cases of ordered 
ground state and disordered one the approach presented 
in this paper can be checked by inelastic neutron scatter- 
ing: corresponding intensity for the singlet-triplet transi- 
tions should have periodicity in the reciprocal space cor- 
responding to the star lattice. This picture is similar 
to observed one in the dimerized spin-Pairls compound 
CuGeOa |23|. In this case inelastic magnetic scattering 
has a periodicity which corresponds to the dimerized lat- 
tice. 

Cases of S > 1/2. — It should be noted that our consid- 
eration of S" = 1/2 KAF can not be expanded directly on 
the cases of larger spins. Although functions presented 
in Fig. 12 where a bold line denotes the singlet state of 
the corresponding two spins, remain eigenfunctions of the 
Hamiltonian for S > 1/2, we have found numerically that 
they are not ground states of the star with S = 1 and 
S = 3/2. All details of calculations are presented in Ap- 
pendix So an other approach should be proposed for 
KAFs with 5 > 1/2. 



tonian of a spin-i magnet in the external magnetic field 
given by Eq. IjlUt . The Hamiltonian parameters are cal- 
culated in the first three orders of perturbation theory 
and are summarized in the Table ^ Within our preci- 
sion KAF has an ordered singlet ground state with all 
the stars in the state given by Eq. JSJ. The ground state 
energy is lower than that calculated in previous finite 
clusters studies. We show that our model can not be 
expanded directly on KAFs with S > 1/2. 

The approach discussed in this paper can be verified 
experimentally on inelastic neutron scattering: the corre- 
sponding intensities for singlet-triplet transitions should 
have periodicity in the reciprocal space corresponding to 
the star lattice. 
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Program " Collective and Quantum Effects in Condensed 
Matter" , the Russian Foundation for Basic Research 
(Grants No. 03-02-17340, 00-15-96814) and Russian State 
Program " Quantum Macrophysics" . 



APPENDIX A: IRREDUCIBLE 
REPRESENTATIONS OF THE GROUP Cev 

The symmetry group Cqv contains six rotations C'^ 
on the angles 27rfc/6 {k = 0,1,. .,5) and six reflections 
which can be written as C*''mi, where ui is operator of 
a reflection. One-dimensional irreducible representations 
can be presented as follows ,2^, ^] : 





-1, 


ui ^ 1; 


(Al) 




- ("1)', 


ui ~ 1; 


(A2) 




- 1, 


ui ^ -1; 


(A3) 




- (-1)', 


Ui ~ —1. 


(A4) 


For two-dimensional ones we have 


MM 





1 

















where two inequivalcnt representations are given hy I = \ 
and / = 2. 



IV. CONCLUSION 

In this paper we present a model for the low-energy 
physics of spin-i kagome Heisenberg antiferromagnets 
(KAFs). The spin lattice can be represented as a set 
of stars which are arranged in a triangular lattice and 
contain 12 spins (see Fig.^. Each star has two degener- 
ate singlet ground states with different symmetry which 
can be described in terms of pseudospin. It is shown 
that interaction between the stars leads to the band of 
singlet excitations which determines the KAF low-energy 
properties. The low-energy dynamics is described Hamil- 



APPENDIX B: STAR WITH 5" = 1 AND 5" = 3/2 

In this Appendix we present details of numerical cal- 
culations shown that functions presented in Fig. |21 where 
a bold line denotes the singlet state of the corresponding 
two spins, are not ground states of the star with 5 = 1 
and S — 3/2 as it was for S — 1/2. 

A simple numerical method for determination of the 
maximum by module eigenvalue of a Hermitian operator 
H (power method ^22,]) was used. It is based on the 
following statement. Let us consider a state of the system 
/ = '^j^Ci'ipi, where the sum may not include all the 



6 



H eigenfunctions. The maximum by module eigenvalue 
Eextr for given / is determined by the expression 



(/|i7»|/> 



(Bl) 



It becomes evident if one notes that = 

The Eq. ljBl|l can be used in numerical calculations 
in the following way. The corresponding expression is 
calculated for n ~ 1,2, ..^nmax- So one can control the 
convergence comparing results with different n. Studying 
a full set of vectors / and taking large enough Umax to 
match the necessary precision one can find the largest by 
module eigenvalue of H . 

In the case of the star the maximum eigenvalue of the 
Hamiltonian is Emax = 185"^ Ji (this energy has the state 
in which all the spins are along the same direction) and 
the energy of singlet states shown in Fig. El is Ess = 



-6S{S + l)Ji. As Emax > \Ess\ for S > 1/2, we have 
to take H = Ho — WI to investigate the lower Hq levels, 
where Tlo is the star Hamiltonian given by Eq. l|T|l. / is 
an identical matrix and W — {E^ax + £'ss)/2 + Ji- Thus 
eigenvalues of H are shifted down according to those of 
Tio by the same value W so as the largest by module 
H eigenvalue to be equal to the TYq ground state energy 
minus W . 



We have not found the ground state energy for the 
star with 5=1 and S — 3/2 by this method because a 
full set of vectors / should be examined for that. This 
operation demands much computer time. But we have 
obtained studying a number of vectors / that there are 
states lower by the energy at least by the value 1.8Ji 
than those discussed above. The method gave Eextr with 
the prescribed precision to the second decimal place at 
nmax = 100 300 depending on / and S. 
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TABLE I: Contributions to parameters of the effective Hamiltonian Eq. IIOII from terms , V"^ and V'^ of the perturbation 
expansion. Interaction J2 has been taken into account in and terms only. Here M is the number of stars in the lattice. 









y3a 

two-stars 


three-stars 


Totals 


J. 





-O.OO7J1 + 0.002 J2 


-O.OlSJi 


O.OlOJi 


-O.OlOJi + 0.002J2 







-O.OOlJi + 0.007 J2 


-O.OOlJi 


O.OOlJi 


-O.OOlJi + 0.007J2 










0.067 Ji 





0.067 Ji 


h 


-0.563J2 


-0.092Ji - O.2I8J2 


-O.I6IJ1 


O.O8OJ1 


-O.I73J1 -O.78IJ2 


AC' 


-O.OO9J2A/' 


-0.768JiA/'+ I.53OJ2A/' 


-O.36IJ1AA 


O.3O4J1A/' 


-0.825JiAA+ 1.521 J2AA 



"This term implies two stars coupling shown in Fig. |3 and three 
stars interaction in the configuration presented in Fig. ISTbt. 
'Correction to the value Co = — 4.5JiA/' for non-interactive stars. 



II 11 




5 5 



FIG. 2: Schematic representation of a star's two singlet 
ground state wave functions (j>i and 02. A bold line denotes 
the singlet state of two neighboring spins, i.e. (IT)J-l) • ~ 

U)JT),)/^/2. 




FIG. 3: a) A cluster where operator of interaction between 
stars has the same symmetry group Cqv as the whole clus- 
ter; b) The only configuration of three stars giving nonzero 
contribution to the third term in the perturbation expansion. 
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FIG. 4: Interactions between two stars: V = Ji(Si S[ + 
S'l^sf ) andV- = J2(S«s(^>+s(i>sf)+sWsf +S«sf ), 
where upper indexes label the stars. The system is symmetric 
under reflection with respect to the dotted line. 



